Phonon effect on two coupled quantum dots at finite temperature 
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The quantum oscillations of population in an asymmetric double quantum dots system coupled 
to a phonon bath are investigated theoretically. It is shown how the environmental temperature has 
effect on the system. 

PACS numbers: 68.65.Hb 

The asymmetric double semiconductor quantum dots system (DQD) is referred as quantum dot molecule due to 
its similar properties to natural molecule. By using self-assembled dot growth technology [lj we can fabricate these 
molecule -like dots, in which the confined electrons can transfer between each quantum dot through tunneling effect. 
Up to now, such system has been studied extensively. Under the influence of an external oscillatory (optical pulse 
or voltage) driving field, one electron can be excited from the valence to the conduction band in one dot, which 
can in turn tunnel to the second dot Recent progress in semiconductor nanotechnology indicate that these 

features have advantages for the application in many quantum devices, such as QD lasers 4j 5J, QD diodesQ as well 
as quantum computing processes [7] [8]. However, QDs are embedded in the surrounding solid matrix, thus the effect of 
electron-phonon interaction during the tunneling is not negligible [9|. Therefore, the environmental temperature will 
have significant influence on the desired QD devices. In this letter, we analyze how the optically driven asymmetrical 
DQD device depends on the environmental temperature. 

The asymmetrical DQD consists of two dots (the left and the right one) with different geometries and has the 
ground state |0) ( the system without excitation), first excited state |1) ( a pair of electron and hole bound in the 
left dot) and second excited state |2) (one hole in the left dot and one electron in the right dot). Any other states 
are effectively decoupled from these three states because the valence band levels of two dots become far off-resonance 
and the hole cannot tunnel to the right dot consequently. Such model can be shown in Fig.??. Considering that the 
single electron is unavoidably scattered by phonons while tunneling between two dots, the Hamiltonian is given by 



H =5>U>01 + n^lOXll + e-^|l>(0|) + T c (|l)(2| + |2)(1|) 



3=0 



+ ~(|1)<1| - |2}(2|) $>(6+ + b k ), 



(1) 



fc=0 



where Ej is the energy of state) j), T c is the electron-tunneling matrix element, tv c is the frequency of the applied field, 
f2(t) = (0\nE(t)\l), where /ds the electric dipole moment, describes the coupling to the radiation field of the excitonic 
transition, E(t) is the optical pulse amplitude. b^(bk) and Uk are the creation (annihilation) operator and energy for 
/cth phonon mode, respectively, g^is the coupling constant determined by the crystal material. 
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FIG. 1: Level configuration of a double QD system. A pulsed laser excites one electron from the valence band that can tunnel 
to other dot with the application of the voltage. We assume that the hole cannot tunnel in the time scale we are considering 
here. 
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Applying a canonical transformation with the generator P [irj| [Tlj 

9k 
2uj k 



5 = (|l)(l|-|2)(2|)^|^(6+-6 fe ), (2) 



we have 



where 



H' = e s He~ s = H' a + H' Il (3) 
H' = e |0)<0| + (e 1 - A)|l)(l| + (e 2 - A)|2)(2| + £ c fc 6+6 fc , (4) 
= O(e^ t |0)(l|e-' 4 + e- i ^ t |l)(0|e' 4 + r c (|l)(2|e 2A + |2>(l|e- 2A ), (5) 



where 



The Hamiltonian in the interaction picture is given by 

if" =e lH « t H' I e- tH '° t = n[e- t( - s - A ^\0)(l\X(t) + e*^- A ^\l)(0\X(t) + ] 
+ T c [e- tUl2t \2)(l\X 2 {t)+e tUl2t \l)(2\X 2 (t) + ], 

where u>ij = e% — £j, S = luiq — uj c is the detuning of the applied field to the the left dot, and 



(7) 



X(t) = exp[- Jr^(K^ kt ~ b k e-^% (8) 

k 

X(t)+=X(t). (9) 

In what follows, we assume that relaxing time of the environment (phonon fields) is so short that the excitons do 
not have time to exchange the energy and information with the environment before the environment returns to its 
equilibrium state. The excitons interact weakly with the environment so that the equilibrium thermal properties of 
the environment are preserved. Therefore it is reasonable to replace the operators X(t), X 2 (t), X(t) + and X 2 (t) + 
with their expectation values over the phonon number states which are determined by a thermal average and write 
the Hamiltonian as 

H eff =Q[e- l(S - A)t \0){l\ +e l(,5 - A)t |l)(0|]e- (A ^ + 5) A 

+ T c [ e -^ 12t |2)(l| +e^ 12t |l)(2|]e- 2(Ar ^ + ^\ ' '" ! 

where A = E^fc(gfc/2wfc) 2 is the Huang-Rhys factor which corresponds to the electron-phonon interactions (l2l|. As 
a result of quantum lattice fluctuations, the exciton-phonon interaction affects our quantum system even at zero 
temperature. Here we have made an assumption that all the phonons have the same frequency, i.e., u>k — ujq, and 
write the phonon populations as N p h — -5^7737 [H I IH 3 • 



Now we proceed to solve the equation of motion for the single electron state vector|0(i)}, i.e., 

ij t m))=H eff m))- (11) 

In general, \4>(t)) is a linear combination of |0), |1) and |2), 

\<j>(t))=c (t)\0)+c 1 (t)\l} + c 2 (t)\2). (12) 
Suppose the initial state of the system is in its ground state |0), i.e., |0) = 1 and |1) = |2) = 0. We have 

ij t Co(t) = ne-^+^c^e-^'^, (13a) 
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FIG. 2: The time evolution of the population in state \2) in the case 5 = A and U12 = 0, for the parameters Q = 0.2ui c , 
T c = 0.1ui c and A = 0.01. Solid line is the result for T = 0, dash line for T = IQuio, dot line for T = 30ojo and dash dot line for 
T = 50lj . 



ijfxit) = fir^^+^dWe'^^ 1 + T c e- 2 ^+5) A C2 (t)e^ 12 *, (13b) 

ij t c 2 {t) = T c e- 2 ^ h+ ^ x c 2 {t)e- tuJl2t . (13c) 

Assuming the optical pulse has a broad square shape, the case 8 = A (the applied optical field is in resonance with 
the zero-phonon line) and u>i2 = (state |1) and |2) are in exactly resonance) yields. 

P (t) = \sin 2 6'cos{®'t) + cos 2 6>'| 2 , (14a) 
Pi(<) = \sin6'sin(G't)\ 2 , (14b) 

P 2 (t) = sin 2 e'cos 2 e l \cos{Q l t) - 1| 2 . (14c) 

where P, = |c,(i)| 2 stands for the possibility of state \i), 9' = ^W 2 + T' 2 , and cos6' = T'J&, with Q' = n e -^ N - h+ ^ x 
and T' c — T c e~ 2 ^ Nph+ ^ x . This analytical result indicates that, in this simplest case, the coherent population oscillation 
are similar to those at zero temperature. But the Rabi oscillation frequency £1 and the electron-tunneling matrix 
element T c are renormalized by the factor of exp[—(N p h + ^)A] and exp[— 2(N p h + |)A] in stead of exp(— ^A) and 
exp{— A), respectively, as shown in FigEJ We note that with the increase of temperature, the Rabi frequency drops 
with the decline of the amplitude. Hence we plot in Fig [3] the average occupation of state \i) [1/too J °° Pi(t)dt] 
(i = 0, 1,2) as function of the temperature, which is obviously show that it is less possible to create the state with 
one electron in the second dot (state |2)) when the temperature is high. 

On the other hand, the environment temperature will also affect the beat pattern, which is caused by the electron- 
phonon coupling, with the rise of the temperature, the beat pattern will decay and the population of the state |2) 
will decrease as well. As the temperature reachs about 50wo, the beat pattern almost disappears. So the environment 
temperature should be considered in the practical applications. 

In conclusion, we have investigated the effect of the environment temperature on an asymmetrical double quantum 
dot driven by an optical pulse. The Rabi oscillation frequency 57 and the electron-tunneling matrix element T c are 
renormalized by the factor of exp[—(N p h + ^)X\ and exp[— 2(N p h + V)X\. With the rise of the environment temperature, 
the population of state \2) declines dramatically and the beat pattern, which is caused by the phonon effect, decays. 
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FIG. 3: Average occupation of state |0) (solid line), state |1) (dash dot line) and state |2) (dot line) as a function of temperature 
in the case 8 = A and W12 = 0, for the parameters Q = 0.2u> c , T c = 0.1tu c and A = 0.01. 




FIG. 4: The time evolution of the population in state |2), for the parameters 8 — U12 = 0, fl = Q.2uo c , T c = Q.1uj c , A = 0.01 
and A = 0.05u; c . (a) is the result for T = 0, (b) for T = 10w , (c) for T = 30w and (d) for T = 50w 
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